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We present a supersymmetric description of the quantum Hall effect (QHE) in graphene. The 
noninteracting system is supersymmetric separately at the so-called K and K' points of the Brillouin 
zone corners. Its essential consequence is that the energy levels and the Landau levels are different 
objects in graphene QHE. Each energy level has a four-fold degeneracy within the noninteracting 
theory. With the Coulomb interaction included, an excitonic gap opens in the zero-energy state, 
while each nonzero energy level splits into two levels since up-spin and down-spin electrons come 
from different Landau levels. We argue the emergence of the plateaux at u = ±(4n — 2) for small 
magnetic field B and at = 0, ±1, ±2n for large B with n natural numbers. 



Introduction: The quantum Hall effect (QHE) is one 
of the most remarkable phenomena discovered in the last 
centuryP,Q|. Electrons, undergoing cyclotron motion in 
magnetic field B, fill Landau levels successively. Each 
filled energy level contributes one conductance quantum 
e^/h to the Hall conductivity a^y The Hall plateau de- 
velops at a^y — iy{e'^/h), where i' is the filling factor. 
It tells us how many energy levels are filled up. Hall 
plateaux have been observed at = 1, 2, 3, ... in the con- 
ventional semiconductor QHE. 

Recent experimental developments have revealed un- 
conventional QHE in graphene0, Q [E The filling 
factors 0,13,01 form a series, v = ±2, ±6, ±10, • • • , where 
the basic height in the Hall conductance step is Ae'^/h 
[Fignja)]. A recent experiment [l(l| has shown a fine 
structure sX v — 0,±1,±4 when larger magnetic field 
is applied. In this Letter, we present a mechanism how 
these unconventional filling factors arise on the basis of 
the supersymmetric (SUSY) quantum mechanics^^l- We 
also elucidate how the graphene QHE is different from the 
semiconductor QHE. 

The low-energy band structure of graphene is de- 
scribed by cones located at two inequivalent Brillouin 
zone corners called the K and K' points. In these cones, 
the two-dimensional energy dispersion relation is linear 
and the dynamics can be treated as 'relativistic' Dirac 
electrons ^3, 0|, in which the Fermi velocity v-p of the 
graphene is substituted for the speed of light. The 
graphene system has the pseudospin degree of freedom, 
where the electron at the K (K') point carries the up 
(down) pseudospin. 

The graphene QH system possesses a unique charac- 
ter that it has a supersymmetry within the noninteract- 
ing theory This follows from the basic fact that the 
intrinsic Zeeman energy is precisely one half of the cy- 
clotron energy for Dirac electrons. It has two important 
consequences [FigQIb)]; the emergence of the zero energy 
state, and the degeneracy of the up-spin and down-spin 
states for each nonzero energy level. Since this holds sep- 
arately at the K and K' points, each energy level has a 
four-fold degeneracy [Fig^b)], and the noninteracting 
system has the SU(4) symmetry. The resulting series 



\s V — ±2, ±6, ±10,- • •. In this paper we focus on the 
Coulomb interactions, and explore how this degeneracy 
is modified. The zero energy state is distinctive, since 
it contains both electrons and holes. Electron-hole pairs 
form an excitonic condensation, producing an excitonic 
gap. Hence, all states become gapful. According to the 
SUSY spectrum, a single energy level contains up-spin 
and down-spin electrons belonging to different Landau 
levels [Fig^b)]. It implies that their wave functions and 
hence their Coulomb energies are different. The Coulomb 
Hamiltonian, projected to a single energy level, possesses 
only the U(l)(g)U(l)(g)Z2 symmetry. Thus, the Coulomb 
interaction leads to the resolution of the four-fold degen- 
eracy into two two-fold degeneracies, each of which has 
the U(l) symmetry. As a result we obtain the new series 
V = 0, ±1, ±4, ±8, • • • . Since all gap energies are pro- 
portional to the Coulomb energy, or ^/B, this series is 
expected to appear for large magnetic field. 
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FIG. 1: (a) The QH conductivity in graphene. The dotted 
black curve shows the sequence v = ±2, ±6, ±10, ■ ■ • , while 
the solid red curve the sequence v = 0, ±1, ±2, ±4, ■ ■ • . (b) 
The energy level and the Landau level. The spin is indicated 
by a solid red (open blue) arrow for electron and hole at the 
K and K' points. This energy spectrum is a manifestation of 
SUSY. The A'^th energy level contains up(down)-spin electrons 
from the (A-l-l)th Landau level and down(up)-spin electrons 
from the Ath Landau level at the K (K') point. 

SUSY spectrum: We start with a concise review of 
the SUSY description of graphene0|. Corresponding 
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to the K and K' points (r = ±), we have two Dirac 
Hamiltonians 

= VF{aa:iPx " tHK^) + Tay{Py - ThKy), (1) 

where Pi = —ihdi + eAi is the covariant momentum, and 



fT, 
(T^ 



(2) 



with cr' the PauH matrix for spins and K = (47r/V3a,0). 
Here, zLK represent the two Brillouin zone corners with a 
the lattice constant. We assume a homogeneous magnetic 
field B = V X A = (0, 0, -B) with B > 0. 

The K dependence is removed from the Hamilto- 
nian by introducing the wave function iprix) = 
ipr{x). Then the Hamiltonian 1^ is expressed as 



rKx~ 



Qd 



(3) 



Q± 

with Q± = Up {cxPx ± ^yPy)- It is diagonalized[lS|, 

= diag. (v/Q^, -VQ^) , (4) 



where the negative component describes holes. 
We consider the quantity 

i7± = Q±g± = vl i~ihV + eAf T ehvla^B, (5) 

where the direction of the magnetic field is effectively 
opposite at the K and K' points. Since this has the same 
form as the Pauli Hamiltonian with the mass m* — 1/ 4u| 
except for the dimension, we call it the Pauli Hamiltonian 
for brevity. The saHent feature of the relativistic Dirac 
Hamiltonian is that its spectrum is mapped from that of 
the nonrelativistic Pauli Hamiltonian. Thus, the energy 
spectrum £„ of the Dirac Hamiltonian is constructed once 
we know the one En of the Pauli Hamiltonian. 

In the Pauli Hamiltonian, the first term is the kinetic 
term while the second term is the Zeeman term. It is fixed 
uniquely as an intrinsic property of the Dirac theory: We 
may call it the intrinsic Zeeman effect. The Landau level 
is created by electrons making cyclotron motion. In the 
conventional QHE, since the Zeeman energy can be con- 
sidered much smaller than the Landau-level separation, 
we may treat it as a perturbation. However, this is not 
the case in graphene. 

The Hamiltonian is a simplest example of the SUSY 
quantum mechanics jllj], where the superalgebra reads 



— 2^Qt, Qt} 



(6) 



with Qr the supercharge (r = ±). The energy eigen- 
values of the Dirac Hamiltonian are found'lJl to be 



7 + T 



and 



F^'' — — F^^ — F^'^ \-hij) \/r 



1 



(8) 



(9) 



for 71 > 0. There exists one zero-energy state only for up- 
spin (down-spin) electrons at the K (K') point [FigQIb)]. 
The existence of the zero energy state is an intriguing 
property of the SUSY theory, where the bosonic and 
fermionic zero-point energies are canceled out^3|. The 
physical reason is that the intrinsic Zeeman splitting is 
exactly as large as the Landau level separation in the 
Pauli Hamiltonian ij^l. The SUSY spectrum tells us that 
it is necessary to make a clear distinction between the 
energy level and the Landau level. This point has been 
overlooked in all previous literatures on graphene QHE. 

The SUSY spectrum dictates that each energy level 
has a four-fold degeneracy [FigQIb)]. The noninteracting 
theory has the SU(4) symmetry. The resulting series is 
j/ = ±2,±6,±10, 

Projected Coulomb Hamiltonian: We include the 
Coulomb interaction to the noninteracting theory, and 
analyze how it affects the four-fold degenerated electron 
states in the Nth energy level. We project the Coulomb 
interaction to the A^th energy level, making a sim ple g en- 
eralization of the lowest-Landau-level projection[lg fa- 
miliar in the conventional QHE. 

For definiteness we expHcitly analyze the TVth energy 
level with iV > 0. The result for TV < is obtained 
by the electron-hole symmetry without any calculation. 
The case = is analyzed in a similar way: See the 
discussion on the exicitonic condensation we give soon 
after. 

According to the SUSY spectrum [FigHtb)], the Nth 
energy level contains both up-spin (down-spin) electrons 
from the (A^+l)th Landau level and down-spin (up-spin) 
electrons from the A^th Landau level at the K (K') point. 
In each energy level there exist four types of electrons 
described by four different field operators 7/)^^(a;). They 
are expanded as 



n 
n 



(10) 



and 



in terms of the wave function 



(7) 



{x) 



AtKx 



{x\N, n) 



(11) 
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and the annihilation operator c^(n) acting on the Fock 
state \N,n) = \N) \n) in the A^th Landau level with n 
the Landau-site index. 

We decompose the electron coordinate x = {x,y) into 
the guiding center X = {X, Y) and the relative coordi- 
nate R = {Rx,Ry), X = X + R, where Rx = —Py/eB 
and Ry = Px/^B with P = {Px,Py) the covariant mo- 
mentum. The projection is to quench the motion in the 
relative coordinate. 

In the A^th energy level, by decomposing the guiding 
center X and the relative coordinate R, the density op- 
erator reads 



arr' 



where b%^, [q) is the projected density fl^ 



(TTT' 

(12) 



D7r' il) 



-E 

27r ^ 



(m|e" 



-i[q+TK-T'K]X 



and F!^^, (q) is the form factor 



(13) 



Fj+ {q) = Fl_{q) = {N\e-^nN), 
Fj+ (q) = fL (g) = {N + l\e-^'i^\N + 1), 
fI_ {q) = Fi+ {q) = {N + lle-^^-J-^^-^lTV), 
f1+ (q) = fI_ {q) = {N\e-''-i+^^^\N + 1). 



(14) 



The set of these form factors is a direct consequence of 
the SUSY spectrum. For instance, F]^_ represents the 
transfer of the up-spin electron {a =t) from the K' point 
(r = — ) to the K point {t — +). In this process an elec- 
tron in the iVth Landau level is moved to the (iV+l)th 
Landau level. Thus, the form factor necessarily mixes the 
A^th Landau level and the (iV+l)th Landau level. This is 
an essential difference from the conventional QHE, where 
only one Landau level is involved. The form factors are 
expHcitly given by using^^ 



{N+Mle'^^lN) = 



ViV! 

^{N + M)\ VV2 
X e 



N 



%q^ 



(15) 



for M > in terms of the associated Laguerre polyno- 
mial. 

The projected Coulomb Hamiltonian is 



rr'a XX' 



(16) 



where 

K^;^;,, (g) = V{q)F^,, {-q)F^'y (q) (17) 
with the Coulomb potential V{q) = e^/47re|<7|. 



The Coulomb Hamiltonian Hpf is invariant under two 
U(l) transformations. 



„T 



(18) 



with two arbitrary constants a and (3. Additionally it is 
invariant under the Z2 transformation. 



(19) 



as corresponds to the fact that the magnetic field is op- 
posite at the K and K' points. The symmetry is quite 
small because of the Landau-level mixing in a single en- 
ergy level, as follows from the SUSY spectrum of the 
noninteracting theory. 
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FIG. 2: A schematic illustration of the energy spectrum, (a) 
Coulomb interactions are neglected. All states are four-fold 
degenerated, (b) Excitonic condensation is taken into account 
at the zero-energy level. The four-fold degeneracy splits into 
two two-fold degeneracies with the gap energy 12811 . Electron 
(solid red line) and hole (dotted blue line) states are sepa- 
rated, (c) All Coulomb interactions are taken into account. 
Each level splits into two subbands. The U(l) symmetry is 
exact in each two-fold degenerated level. 

Excitonic condensation: It is necessary to treat the 
zeroth energy level {N = 0) separately from the others, 
since it contains both electrons and holes. They come 
from the A^th Landau levels with A'' = ±1. Due to the 
Coulomb attraction, electron-hole pairs are expected to 
make bound states and condense into the spin-singlet ex- 
citonic states. 

It is sufficient to investigate only electron-hole pairs 
at the K and K' points separately, because an exciton 
composed of an electron and a hole belonging to different 
Dirac cones is fragile. At the K point their field operators 
are 



n 



(20) 



The dominant electron-hole interaction is 



Hn^ttJ d^qViq)i;li-q)M-Q)^liq)MQ)- (21) 
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The analysis of this Hamiltonian is similarly done as in 
the BCS theory. However, the kinetic term is absent 
since it is quenched in each Landau level. This simplifies 
the analysis considerably. In the mean-field approxima- 
tion, the singlet excitonic gap function satisfies the gap 
equation, 

A(fe)= /rf2^V(fc'-fc)e-'«('='-'=)'/2tanh^^, (22) 
J "ikBT 

where is the Boltzmann factor, and the suppression 
factor e^'sC" /2 jj^s arisen from the projection of the 
Coulomb interaction. In the limit T ^ 0, the zero- 
momentum gap A (0) is given by 

A(0)|j,^o = 7rV2^(eV47re/B) ■ (23) 

The excitonic condensation resolves the electron-hole de- 
generacy in the zero-energy state. As a result the four- 
fold degenerated levels split into two two-fold degener- 
ated levels [FigEIb)] with the gap energy Ij23|l . This leads 
to a new plateau at i/ = 0. 

Pseudospin asymmetry: We can treat electrons and 
holes separately since a gap has opened between the elec- 
tron and hole bands [FigEtb)]. Each energy level is four- 
fold degenerated within the noninteracting model with 
the SU(4) symmetry. However, the Coulomb interaction 
breaks it explicitly into U(l)(g)U(l)®Z2. We now argue 
that plateaux emerges at = ±l,±2n by this symmetry 
reduction [Figl^c)]. 

We take the Hartree-Fock trial function for the A^th 
energy level {N ^ 0) by requiring the symmetry, 

l*]v> - n ("^^+ + (")) |0> ' (24a) 

n 

|$jv> = n ("^^+ (") + ^^'^-^ ' (24b) 

n 

with -I- jw'^p = 1. They transforms properly under 
the U(l) transformation Ijl8|l . and exchange themselves 
under the Z2 transformation Ijlijll . 

It is easy to determine u'^ and v"' by minimizing the 
Coulomb energy {^'J^\Hc\^n) ■ The result says 

= v'^ — e*^ sin a, = ^ e^'^ cos a, (25) 

where a is a certain constant given in terms of integrals 
over various form factors together with the Coulomb po- 
tential. The arbitrary phase assures the U(l) symme- 
try. It follows that Im^^I > \v^, leading to the pseudospin 
asymmetry, because the Coulomb energy of an electron 
in higher Landau level is lower. These two states are 
degenerate, 

($]^|i7c|$]v> = (*]vl^^cl$]v). (26) 
due to the Z2 invariance Ijl9|l . 



To study the zeroth energy level, we take the trial func- 
tion 

1*0) = n ("c? (n) + vc}J (n)) |0) (27) 

n 

with jitp -I- |wp — 1. By minimizing the energy we find 
either u — or v ^ 0, refiecting the Z2 symmetry. 

We have so far taken the intrinsic Zeeman effect into 
account. However, there may be an additional Zeeman 
effect in graphene, which make these two states split ex- 
plicitly. Even without such an extrinsic Zeeman effect, 
driven by the Coulomb exchange interaction, the spon- 
taneous breakdown of the Z2 symmetry turns the system 
into a QH ferromagnet (iflj . In any case the excitation gap 
is of the order of the typical Coulomb energy. This ex- 
plains the emergence of plateaux at v = ±1, ±2, ±4, • • • . 

Discussions: We have presented a SUSY description 
of the QHE in graphene. It has an important conse- 
quence that the energy levels and the Landau levels are 
quite different objects in graphene [FigQIb)]. One en- 
ergy level contains electrons coming from two neighbor- 
ing Landau levels. This is evidenced in the Coulomb 
energy through the form factor ifTiji . We have derived 
the Coulomb Hamiltonian ((Tfill . as dictated by the SUSY 
spectrum. 

Our formalism is considerably different from those as- 
sumed by previous authors. Let us comment on them. 
First of all, Alicia and Fisher ji^l have not taken into ac- 
count the projection of the Coulomb interaction. They 
have also ignored the Landau-level mixing in a single en- 
ergy level. Nomura and MacDonald^J] made only an 
explicit analysis of the SU(4) symmetric term, which is 
extracted from lfT2|l as 

pT'\q) = llFNiq) + FA.+i(g)} ^ (q) , (28) 

ar 

where we have set FNiq) = {N\e-'i^\N). We should 
note that the SU(4) noninvariant term is as large as the 
invariant term. This is essentially different from the con- 
ventional QH system, where the noninvariant term can 
be made arbitrarily small by controlling the external pa- 
rameters such as the layer separation d, the tunneling gap 
AsAS and the magnetic g-factor. Finally, the spin degree 
of freedom has been ignored by Goerbig et al.[23, where 
F/v (g) and Fjv+i (q) appear only in the SU(4) symmetric 
combination in the density operator. 

We have started with the SU(4) symmetry in the 
noninteracting theory, where gaps open at i' = 
±2, ±6, ±10, • • • . The Coulomb interaction breaks it ex- 
plicitly into the U(1)(E'U(1)(8)Z2 symmetry. As a result, 
new gaps open at v — ±1, ±4, ±8, • • • . The remaining 
problem is whether a further resolution of the degener- 
acy may occur. 

Here, we recapture a similar problem in the spin- 
frozen bilayer QH systemQ- Let us ignore the tunneling 
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gap AsAS- Then, there is the pseudospin SU(2) sym- 
metry in the noninteracting theory, which is broken to 
the U(l) symmetry explicitly by the capacitance effect. 
Here, the exchange Coulomb interaction generates the 
pseudospin wave, which is the Goldstone mode associ- 
ated with spontaneous breakdown (SSB) of a continuous 
symmetry, and turns the bilayer system into the pseu- 
dospin QH ferromagnet[^^. Based on this analogy No- 
mura and MacDonald!2lj has concluded a spontaneous 
development of the SU(4) QH ferromagnet in graphene. 

We question why the SSB of a continuous symmetry 
is possible at finite temperature in spite of the Mermin- 
Wagner theorem in the conventional QHE. This is be- 
cause there exists the tunnehng gap Asas 7^ in the 
actual bilayer system, which gives a gap to the Gold- 
stone mode. This is not the case in the graphene QHE. 
Graphene is an ideal two-dimensional system, and fur- 
thermore the U(l) symmetry is exact. There is no ex- 
ternal parameter which breaks it explicitly to make the 
Mermin- Wagner theorem inapplicable. We conclude that 
the SSB of the U(l) symmetry cannot occur in graphene, 
so that the plateaux at = ±3,±5,±7, ••• will not 
emerge. 

The author is grateful to Professors T. Ando, Y. lye 
and Z.F. Ezawa for fruitful discussions on the subject. 
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